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We show that the Mellin summation technique (MST) is a well defined and useful tool to 
compute loop integrals at finite temperature in the imaginary-time formulation of thermal field 
theory, especially when interested in the infrared limit of such integrals. The method makes use 
of the Feynman parametrization which has been claimed to have problems when the analytical 
continuation from discrete to arbitrary complex values of the Matsubara frequency is performed. 
We show that without the use of the MST, such problems are not intrinsic to the Feynman 
parametrization but instead, they arise as a result of (a) not implementing the periodicity 
brought about by the possible values taken by the discrete Matsubara frequencies before the 
analytical continuation is made and (6) to the changing of the original domain of the Feynman 
parameter integration, which seemingly simplifies the expression but in practice introduces a 
spurious endpoint singularity. Using the MST, there are no problems related to the implementation 
of the periodicity but instead, care has to be taken when the sum of denominators of the 
original amplitude vanishes. We apply the method to the computation of loop integrals appearing 
when the effects of external weak magnetic fields on the propagation of scalar particles is considered. 
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I. INTRODUCTION 

At finite temperature, unlike in vacuum, momentum 
dependent loop integrals in general depend separately on 
the time (po) ^^nd space (p) components of the momen- 
tum = {pq, p), since Lorentz invariance is lost due to 
the presence of the medium. As a consequence, the lim- 
iting behavior of these integrals, as the momentum com- 
ponents approach given values, may depend on the way 
the limit is taken. For instance, in one loop self-energy 
calculations, the infrared limit {pq ~ 0,p = |p| — > 0), 
that accounts for the plasma screening properties, does 
not necessarily coincide with the static limit {po —^0,p — 
IpI = 0), that accounts for the long wavelength plasma 
oscillations. In physical terms, this non-analyticity as 
^ is due to the cut structure of the self-energy at 
finite temperature, where branch cuts appear represent- 
ing scattering processes not allowed in vacuum. 

The above behavior was originally not fully recognized 
since there were results indicating that the aforemen- 
tioned two limits commuted when cornputing the real 
part of the self-energy in a (j)^ theory [U, Q • In the imag- 
inary time formalism (ITF), the problem with these cal- 
culations was traced back to an incorrect analytic con- 
tinuation when the discrete frequency takes on arbitrary 
complex values [3]. Since the erroneous result was ob- 
tained from an analysis based on the use of the Feynman 
parametrization at finite temperature, it is often thought 
that such parametrization is also endemic to the source 



of the error. 

In the ITF, the introduction of the Feynman 
parametrization for the computation of loop integrals 
containing two propagators can be avoided since there 
are techniques that allow to perform the sum over Mat- 
subara frequencies in a straightforward manner. How- 
ever, it has been recently shown that when studying the 
influence of weak magnetic fields over physical processes 
at finite temperature, the loop integrals that appear in- 
volve products of powers of two or more propagator-like 
denominators Although it is possible to generalize 
the standard techniques to carry out the sum over Mat- 
subara frequencies from the product of two propagators 
to the case a of product of powers of these, the calcu- 
lations become extremely cumbersome. It is therefore 
desirable to have a more direct method to perform these 
calculations. One of such is the Mellin summation tech- 
nique (MST). The method [1] calls for the use of Feyn- 
man parametrization, which allows to condense products 
of powers of propagators into a single propagator-like fac- 
tor raised to some power. This is particularly useful when 
one seeks an answer in terms of a power series involving 
a small parameter, for instance, at a high temperature 
T, the ratio m/T, where m is the particle's mass. 

In light of the well known mishaps with the use of the 
Feynman parametrization [sl] in finite temperature cal- 
culations, it is important to establish that the MST in 
the ITF has the correct analytical properties when the 
discrete Matsubara frequency is continued to arbitrary 
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complex values. In this work we undertake such study. 
We perform an explicit calculation of the one-loop self- 
energy in a (jr" theory as a guiding tool to find out how 
the Feynman parametrization should be used in the ITF 
of thermal field theory. The work is organized as fol- 
lows: In Sec. In] we give a brief summary of the previous 
analyses that have dealt with this problem and compare 
their results. In Sec. IIIIl we give a detailed derivation of 
the standard calculation using conventional techniques. 
We explore the case when the external momentum ap- 
proaches zero and in particular, give the explicit result in 
the infrared limit. As we want to compare with the MST, 
the answer is given in terms of an expansion in powers of 
m/T for the case when T ^ m. We also point out the im- 
portance of considering that the external frequency takes 
on discrete values and therefore implementing the period- 
icity of the resulting expressions before taking the analyt- 
ical continuation to arbitrary complex values. In Sec. IIVI 
we perform the calculation using the MST which involves 
the use of the Feynman parametrization. We underline 
the importance of carrying out the integral over the Feyn- 
man parameter x € [0, 1] to avoid the appearance of spu- 
rious endpoint singularities. We carefully show how the 
simple Feynman parametrization has to be corrected to 
account for the case where the sum of denominators in 
the Feynman formula vanishes and emphasize that this 
correction term accounts for the whole dependence on the 
way the momentum approaches zero, in agreement with 
the analysis in Ref. p']. In Sec. |V]we apply the results 
for the computation of integrals describing the self-energy 
of a neutral scalar interacting with charged ones in the 
presence of a magnetic field. We finally present our con- 
clusions and give an outlook in Sec. IVII We leave for the 
appendices the demonstration of important intermedi- 
ate results and alternative derivations of the calculations 
arising in the discussion of Sees. IIIIl and IIVI 



II. NON-ANALYTICITY AND MISHAPS WITH 
FEYNMAN PARAMETRIZATION 

The problem of the non-analyticity of thermal field the- 
ory calculations as the momentum components approach 
zero, has been analyzed by several authors. The land- 
scape of findings is, at first sight, rather blurred since 
there are many details in the calculations that are sources 
of more extended discussions. Among these we can men- 
tion: the implementation of derivative expansion tech- 
niques; the validity of perturbative and derivative ex- 
pansions exchange and the implementation of the exter- 
nal bosonic field periodicity. Other studies are concerned 
with the correct analytic continuation to arbitrary com- 
plex values of the external frequency; the soundness of 
some redefinition of variables inside potentially divergent 
integrals; the physical interpretation of the imaginary 
part of the thermal bubble and the use of the Feynman 
parametrization. 

In this work, our main purpose is to show that the 



use of Feynman parametrization, within the MST, is a 
well defined procedure. Nevertheless, it is worth paus- 
ing to summarize, from a wider perspective, what has 
been found in the context of the non-analyticity of ther- 
mal self-energies at the origin in calculations that do not 
resort to the use of the MST. 

In general, this problem has been dealt in terms of 
perturbative and non-perturbative approaches and both 
in the ITF and the real time formulation (RTF) of ther- 
mal field theory. The coincidence between perturbative 
and non-perturbative calculations is sometimes taken as 
a guide to decide on the correctness of the approach. 

The discrepancy between the results in the infrared 
and the static limits caused a great deal of confusion, 
prompting a number of possible explanations. These 
ranged from assigning validity only to the infrared limit 
as a genuine result in thermal equilibrium JJ , passing by 
suggesting that it is not necessary to assume that the 
external field is in thermal equilibrium [H, d, , to dis- 
missing the non-analyticity by claiming that this is not 
present in an exact solution to the slow motion approxi- 
mation of the Green's function P, 0] ■ In Ref- [H , the cal- 
culation in the ITF is done by extending the external fre- 
quency to the whole imaginary axis and then analytically 
continuing it to the entire complex plane. With this pro- 
cedure, the periodicity of the functions in the external fre- 
quency -that was present before analytical continuation- 
is lost. This is how the erroneous result, that the infrared 
and static limits coincide, is obtained. This result seem- 
ingly confirmed the one in Ref. which was performed 
in the RTF. Reference @ points out that truncating the 
derivative expansion at the beginning of the calculation, 
either by keeping only the constant term [To| or at higher 
order [l|, 0] , gives misleading results, in the first case, be- 
cause the operator nature of the background field is lost; 
in the latter, because the periodicity is not considered. 

However, Weldon Q showed that the results in Refs. [l|, 
^ go wrong, performing the calculation both in the ITF 
and the RTF. In addition to providing physical argu- 
ments for the inequivalence of the infrared and static 
limits he demonstrated that in the RTF calculation, the 
use of the Feynman parametrization, as is commonly im- 
plemented in T = calculations, needs to be corrected. 
The correction accounts for the fact that the real time 
Feynman amplitude is not the boundary value of a single 
analytic function and thus it is necessary to perform one 
calculation for the real and another one for the imagi- 
nary part. However, for the discussion concerning the 
ITF, Weldon argued that starting from an expression 
that uses the Feynman parametrization, the analytical 
continuation is not unique and leads to a function con- 
taining branch points and an endpoint singularity that 
need to be removed by the addition of an extra term. 
The argument is based on an expression where the inte- 
gration interval for the Feynman parameter x has been 
changed from x G [0,1] to x € [0,1/2] which uses the 
symmetry of the integrand about x = 1/2, before the 
analytical continuation is implemented. 
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In what follows, we will show that the Feynman 
parametrization can be implemented in the ITF (within 
the MST) without introducing spurious branch points 
and endpoint singularities. The key ingredients are the 
implementation of the periodicity in the expressions be- 
fore the analytical continuation, the use of the original 
integration interval for the Feynman parameter x and the 
accounting of the extra term that corrects the original 
Feynman formula when the sum of denominators van- 
ishes and which happens naturally in the MST. Before 
proceeding with this analysis, it is convenient to set the 
stage and perform an explicit calculation using the stan- 
dard technique in the ITF to have a reference to compare 
with the result obtained after introducing the MST. 



III. STANDARD CALCULATION 

We start the discussion with the explicit expression 
for the one-loop self-energy of a scalar field (p with a self- 
interaction of the form Xff)^ . This is given by 
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The standard calculation 's'l is done by first summing 
over the Matsubara frequencies. This is most easily ac- 
complished by using Eq. ^ to write Eq. ([T]) as 
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The sum over the Matsubara frequencies can be com- 
puted by means of the identity 
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which yields the expression for the self-energy 
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rcothf^V.coth^^)], (5) 



where we used that coth(a; 4- Ht^) — coth(a;), in ac- 
count of the fact that p^i = 2iTTlT and furthermore that 
scoth(sx) = coth(a;). We emphasize that this is an im- 
portant step necessary to implement the periodicity in 
the expression and that eventually allows the analytic 
continuation of the result to arbitrary complex values 
of poi. As discussed in Ref. @ when this condition is 
not taken, the result cannot be interpreted on physical 
grounds. To stress the importance of this point, we show 
in Appendix I that, had this condition not been taken, 
the eventual analytical continuation would have lead to 
an erroneous result that needs to be corrected precisely 
by the addition of the function Ilg found in Ref. [3| . Sim- 
plifying Eq. (O we get 
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Note that upon the change of variable k — p ^ k, the sec- 
ond term in Eq. ([6]) reduces to the first one and thus the 
complete expression for the self-energy is twice the first 
term in the above equation. Carrying out the angular 
integration we get 
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At this point we take the analytical continuation in poi 
from discrete imaginary values to arbitrary complex ones, 
Poi ^ Po and explore the limiting behavior of Eq. ([7]) as 
the momentum components of the vector = (po; , p) 
approach zero. We specialize to the case where po is real. 
Since the result depends on the way the limit is taken, 
we first set pq = ap 
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and take the limit p ^ by expanding the logarithm 
around p = 0, 
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After carrying out the sum in Eq. Q we get in the limit 
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By using the identity 



coth( ^ ) = il + 2n{Ek)), 



(11) 



where n(Ek) is the Bose-Einstein distribution, we can 
separate the vacuum and thermal contributions of the 
above equation. Keeping only the thermal part we obtain 
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We now follow Ref. to find the explicit expression for 
n in the infrared limit at high temperature. The result 
for arbitrary a is given in Appendix II. Setting a = in 
Eq. dm) we get 
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where in order to write the integral in d— dimensions, we 
have first extended the integration domain from [0, oo] to 
[—00,00] and thus multiplied by 1/2. The extension of 
the integral in Eq. ([T^ to d-dimensions represents a way 
of handling the infinities involved in the explicit compu- 
tation and we should keep in mind that in order to make 
contact with Eq. (|12p . the limit d ^ \ will be eventually 
taken. The extension to d— dimensions also calls for the 
introduction of the mass scale /i. The angular integration 
in Eq. can be done straightforward and the result is 
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where /? = 1/T, Eq. IM]) can be written as 
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Upon the change of variable 
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The integrals in the last expression are well known and 
can be expressed in terms of ratios of gamma functions 
r, namely 

\zil-zrz^^^-if^plh±21. (19) 
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Separating the term with n = in the sum and keeping 
in mind that the terms in the sum are even powers of n 
we get 
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The first and third terms within the curly brackets in the 
right hand-side of the above equation have a singularity 
when d = I that should be isolated. The singularity in 
the first term arises as the argument of one of the gamma 
functions vanishes. The singularity in the third term is 
less obvious and we concentrate on it. Defining 
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where C is the Riemann zeta function, which has a simple 
pole at d = 1. In the high temperature limit T ^ to, we 
can approximate the above expression as 

(27rT)''-2C(2 - d) 
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Substituting Eq. 

n^(o,p^o) = 



into Eq. ([2T|) we get 
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We now set d = 1 — 2e and make a series for e 0. 
The e-poles cancel and the expression for if^ at high 
temperature and in the infrared limit is 
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IV. MELLIN SUMMATION TECHNIQUE AND 
FEYNMAN PARAMETRIZATION 



tions in the ITF. The technique resorts to applying a 
Mellin transform over the discrete frequency to the ex- 
pression involving the sum and afterwords applying the 
inverse transform to obtain an identity. In this fashion, 
the calculation of the sum becomes the easiest part and 
the problem reduces to computing the McUin transform 
and its inverse of the remaining expression. 

We now show how the combined use of the MST and 
the Feynman parametrization leads to the same result as 
the standard calculation obtained in Eq. (pS)) . For this 
purpose, we start from the expression for the self-energy 
in Eq. ([T]) separating the sum over Matsubara frequencies 
as 



+ |/| 
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This expression has the advantage of separating the sum 
into pieces where the frequencies involved have a definite 
sign from the one where the frequencies have a mixed 
sign. The former is suited for the application of the MST 
since this last is an integral transform over a continuous 
variable restricted to the positive real axis [see Eqs. ([5^ ]. 
Since the combination of frequencies appears as a square, 
all that matters is that this has a definite sign, either 
positive or negative. For the latter this is not possible. 
Nevertheless the calculation can be performed making 
use of the Feynman parametrization. In addition, notice 
that by transforming the original discrete frequencies into 
a continuous variable, there are no problems associated to 
the implementation of periodicity conditions. However, 
care has to be taken when the sum of denominators of 
the original Feynman amplitude vanishes, leading to a 
correction term, as discussed by Weldon in Ref. In 
what follows we analyze these contributions separately. 



A. Definite sign frequencies 

To begin, let us concentrate on the first two terms aris- 
ing from the separation of the sum over Matsubara fre- 
quencies in Eq. (p6| and define 
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The MST is a useful tool to compute infinite sums [Tl| In Eq. D 
of the form encountered in Finite Temperature calcula- 



where we have introduced the Feynman parametrization 
and thus the integral over the Feynman parameter x. 
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with LOn = 2TTnT and oji = 2-kIT. We shift the three 
momentum integration variable k — > k — (1 — a;)p and, 
after making the appropriate renaming of the summation 
index, the expression for Hi can be written as 



)2 el 



(28) 



where S± are defined as 



n=0 
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with 
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and 
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The quantities uj± wiU become the variables over which 
the Mellin transform is computed, that is to say u!± y 
in Eq. ([5^ . Notice that the expression for / involves 
an integral that for later purposes has been extended to 
d-dimensions, and thus the need to introduce the mass 
scale ^. In order to make contact with Eq. ([77|) . this 
time, as opposed to the discussion after Eq. (fT^ . we will 
be interested in taking the limit d ^ 3. 

We perform the sums S± by means of the Mellin sum- 
mation technique. In general, the Mellin transform pair 
fiy), M[f;s] is given by 



M[f;s]=j;^y^-^f{y)dy, 

fiy) = ^j:-:::y-'M[f;s]d. 



a < Re(s) < (3 
a < c< P, 



(32) 



where a and (3 are determined by the condition that the 
first of the integrals in the above equations converges at 
y — and y = oo, respectively. The variable y contains 
all the dependence on the summation variable n and is 
treated as a continuous variable. By expressing /(y) as 
an inverse Mellin transform, we can then perform the 
summation over n. The problem reduces then to find- 
ing the Mellin transform and its inverse of the remaining 
expression. 

From Eq. ((30)) . it is easy to see that a — Q and /3 — 
4 — d. In terms of their Mellin transforms, S± can be 
expressed as 
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The sum over n can be explicitly evaluated, yielding 
" 1 
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where C('^j^) is the modified Riemann zeta function. 
To find out the Mellin transform of /, we observe that 
the integrand of the first of Eqs. ([32| . with / given by 
Eq. (|30p . can be thought of overall as an integral in 
(s+(i)-dimensions of a function of the square of an (s-l-ii)- 
dimensional vector 



dim d-dim 
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Such integrals are well known [ij] and the result, after 
compensating for the volume of the solid angle when ex- 
tending the integral from d to (s ~t- (i)-dimensions, is 
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where F is the gamma function. Combining the results 
in Eqs. (p4|) and p6|) the explicit expression for S± is 
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In order to perform the integral over s in Eq. (j37p . we 
notice that it is necessary to know whether the term 
[m^ + x{l — x){ujf + p^)]/(27rT)^ is larger or smaller 
than one. For the present purposes where we work 
in the high temperature limit and want to explore the 
analytic properties of 11 near the origin, we see that 
[m^ + x{l - x){iof +p^)\/{2iTTf < 1. Notice that this 
assumption limits the range of values of the external in- 
dex / to be / = 0, ±1. Taking c? = 3 — 2e, e — > 0+, we can 
choose c such that 1 < c < l-f 2e, in order to both, com- 
ply with the upper bound requirement for the existence 
of the Mellin transform, Eq. ([5^ . and to avoid the pole 
of C at s — 1. Therefore, the integration contour can be 
closed to the right by a half-circle at infinity. The only 
singularities within the integration contour are those of 
F[2 - (3 - 2e + s)/2], namely when s = 1 + 2e + 2k, 
k = 0,1,2, .. . and the integral over s in Eq. ([57]) can be 
computed by means of the residue theorem, yielding 
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Notice that for fc > 0, C(l + 2fc + 2e, |l| + 1 ± (1 - x)l) has 
no singularities since |Z| + 1 ± (1 — x)l 7^ 0, —1, —2 . . . For 
fc = 0, the singularity is regulated by e. Also, the factor 
[m^ + x(l — x){u!f + p'^)]'' is non-singular. Therefore we 
can analytically continue the functions S± from discrete 
to arbitrary complex values iui — 2'niTl — > po- Upon 
this analytic continuation 



■ Po 

—I 

2ttT 



2ttT 



(39) 



We now explore the limiting behavior of Eqs. ((38)) as 
the momentum components of the vector ~ (poiP) 
approach zero. Since the result depends on the way the 
limit is taken, again we restrict ourselves to real po values, 
set pq — ap and take the limit p for the argument 
of the C function. 

C[l + 2(fc + e), 1 + {\a\ T «(1 - x)a)p/2TrT] 
C[l + 2(fc + e)]- 

(|a| T *(1 - x)a){l + 2{k + e))C[2 + 2{k + e)] (^) . 

(40) 

The integration over the Feynman parameter involves the 
computation of the integral 



Z± 



dx 



{C[l + 2{k + e), 1 + (|a| T *(1 - x)a)^] 
' x{\ — x){a^ 



P 



2ttT' 



l)p' 



(27rr)2 



(41) 



As we are interested in evaluating the result near the 
origin and in the high temperature limit, we make use of 
the expansion of ( in Eq. into Eq. ([1T|) and evaluate 
for the first two terms in the series expansion of Eq. ([38)) . 
namely k = 0, 1. Defining 



Z± 



k=0 



(42) 



we get 



= C[l + 2e] 

- (|a|T*a/2)(l + 26)C[2 + 26](^) 

+ 0(/) 

Zi = C[3 + 2e] ( 



V27rT. 

(|a| T«a/2)(3 + 2e)C[4 + 2e] 
Oip'). 



27rT/ \2nT 



—) 
InTJ 



(43) 



Notice that in the limit p ^ 0, the result coming from the 
sums S± is independent of a. Combining Eqs. 



and ([43]) into Eq. (128]), we obtain 

ni(o,P-o) = ^ 

2e ^V2\/^T 



7B m2C(3) 
2 



87r2r2 



(44) 



B. Mixed sign frequencies 



We now turn to the computation of the third term 
arising from the separation of the sum over Matsubara 
frequencies in Eq. (pS)) . We define 



^2{P01,P) 



A2 
2 



d^k 

(2^ 



dx 



xujiY + yf 

2\ 



(45) 



where y — (k ~ xp)^ — a;(l — a;)(Po/ ^ P'^) + 
order to perform the sum over Matsubara frequencies, it 
is convenient to note that the degree in the denominator 
of Eq. ([45]) can be reduced since 



1 



d 



1 



XUJi 



dm? [(cj,, 



d 



XLOl) 

is 



-y] 



5m2 ^ 2wi/2 

s=±l " 

^ U 1/21 ' ("^^^ 

[(Wn - xuji)+ isy'-'^l 

where in the last line we have resorted to partial frac- 
tioning. With this reduction, the sum over Matsubara 
frequencies can be easily performed and we get 

'YJ (2^ 



n2(po;,p) 



dx 



d 



1 



dm? ^—^ Av^ 

s=±l " 



/2 



IS 



IT 

is 



IS 



TT 



coth 



sxpoi 



TT 



2T 

sxpoi 
2T 



sxpoi 




2T 


2T 


yl/2 1 




2T 


) 




)} 


2T 





(47) 

where 'f/' is the digamma function. Before proceeding 
further, it is convenient to note that given a function F 
with argument {spt^x + y{x, m)^^"^) /2T , the identity 



d 



F 



( SBh 
\ 21 



l£ 

2T 



2T 



y{x, to)^/2 



d 


p ( spoix 
^ \ 2T 


~^ 2T J 


dx 


2y{x,m) (^spoi 4 


1 dy(x,m)\ 




2iy(a;,m)i/2 dx Jj 



(48) 
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is satisfied. Using Eq. (|48|) into Eq. (|47)l . the integration 
over the Feynman parameter becomes trivial and we get, 
after performing the angular integration 



Uiipohp) = - 



E 



coth 



IS , 

+ — ?A 

TT 

is 

- 

TT 

is , 

+ 

TT 

is 

- —ijj 

TT 



In 

Ek 
2T 



pIi 



- 2spQiEk + 2kp 
Pqi ~P^ ~ '^spoiEk - 2kp 



coth 



2T 



spoi 
2T 



is 


'Ek 


spoi 


) 


TT 


2T 


2T 




is 


'Ek 


spoi 




TT 


2T 


2T 


) 



is Ek 
1+ + 



(49) 



Notice that all along this part of the calculation, we have 
distinguished between I and poi which are in principle 
related through po/ = 2iTTlT. The reason is that the 
first term arises as a consequence of our treating this 
partial sum over Matsubara frequencies as having \l\ in 
the limits of the summation index, whereas the second 
is the discrete value taken by the external energy. If we 
now use that p^i = 2iTTlT we notice in particular that 
the periodicity in Eq. (|49p , coming from the argument of 
the second coth, gets accounted for. However, there is 
an extra term that exhibits periodicity in pm and that is 
not that evident from the expression in Eq. (HH). This 
comes from the difference of the functions '0 that have 
Pqi in their argument. This periodicity can be evidenced 
by resorting to the identity 



'ij}{x + iy) — il'{x — iy) — 



2iy 



k=0 



y2 + + ky 



i TT coth[7ry] 



2iy 



k=l ^ 



(50) 



After this simplification which allows to account for the 
periodicity in pQi , the equation can be analytically contin- 
ued from discrete to arbitrary complex values oipoi ~* po, 
since Eq. (|49|) is free from singularities. As before, we 
explore the limiting behavior of Eq. (|49p as the momen- 
tum components of the vector = {po,p) approach 
zero. Again, we restrict ourselves to real values of po, set 
Po — ap and take the limit p —^ and we get 



A2 



dk 



2{2tt)' 





TT 

2m 



1 



(51) 



which is independent of a. Upon combining the results 
of Eqs. (Hi]) and ^5^, we obtain 



n(0,p^O) 



A2 



4(27r)2 



— 1 ] f 

m ^ 2e^ \2^T 



IE m\(3) 
2 87r2T2 



(52) 



C. a dependence 

We now proceed to discuss the a dependence of the 
result since, as we have seen in Sees. IIVXI and lIVBi this 
dependence is absent in the terms calculated so far. 

As it is shown in Ref. 0, the usual Feynman 
parametrization formula at finite temperature has to be 
corrected when the sum of denominators can vanish. The 
correct expression is in this case 



-D1-D2 



dx 



Atti 



[(1 - X)D2 - 
6{Di + D2) 



xDi] 



(53) 



where the first term is taken as the principal value. Since 
within the MST, the Matsubara frequencies in the sum 
are first treated as continuous upon the Mellin transform, 
and an analytic continuation is required, there is the pos- 
sibility that the second term contributes. We proceed to 
show that this is indeed the case and that this last term 
carries the full a dependence of the result. In Appendix 
II we show that the term here computed coincides with 
the one obtained by using the standard procedure. 

We first look at the contribution from the positive sign 
frequencies which can be written as 



n 



Ji>0 



47riT^ 

n=l 



d^k 



(27r)3 

^ S[ic.,,-u;,r + El_^+u;l + El] 
{uj,,-u:iY + El_^-u:l-El 

Upon the change of variable 

k + p/2. 



(54) 



(55) 



and the introduction of the Mellin transform and its in- 
verse, Eq. ([54|) can be written as 



V ^ / ^Jc-ioo Jo 

d^k 5[{-iu - ojiy - h 



duu^ 



{2ttY 



{—iu — LOiY + — 2p • k 



(56) 



where the Mellin transform has been taken from the dis- 
crete Matsubara frequency a;„ — > — iu, in anticipation for 
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the taking of the analytical continuation iuji — > pq. Also, 
for convergence of the integral, < c < 3. Taking the 
analytic continuation iuji po, setting po — ap, the in- 
tegral over k can be performed straightforward. In the 
limit p — *■ we get 



(Pk S[—{u — ap)"^ 



(2^)2 



In 



~(u — ap)"^ + _ 2p • k 
\/ — + au 
— \Jv?-— vn? + au 



9 ? 

amp 



(57) 



It is easy to show that the potentially dangerous first 
term in the above equation for p — > is canceled from a 
similar contribution arising from the sum over negative 
Matsubara frequencies. We thus just concentrate in the 
second term of Eq. (|57p . We give explicit results for the 
case a > 1; the case a < 1 can be worked out by resorting 
to the transformation formulas for the hypergeometric 
function. Here we jut point out that when < a < 1, 
the high temperature expansion contains an imaginary 
part. Integration over the variable u gives 



2 (1 - a2 



3-s 4-s 1 



2 ' 2 ' 1 - a2 



r(^j 2 ' 2 ' i-a2 

r(§)r(^) (-1)^/2 



3-s 4-s 1 



+ i 



X 2F1 



r(i) (l-a2)s/2-l 

sis 1 
2' 2' 2' 



(58) 



where 2F1 is a hypergeometric function. The remaining 
s-dependent factor comes from the summation over the 
Matsubara frequencies in Eq. ((56|) which yields 



Therefore, the integral over s involves the terms 

rri 2 />C+ZOO 

1 Q -^TT ' 



Li = 



m 2{2m) 

, 3-s 4 
2F1 1, 



ds 



2ttT J (l-a2)r(^) 
1 



2 ' 1 



c-i-ioo 



ds 



*m 2(2^1) 

r(f)r(V) 



-imV (-l)'^/2C(s) 

2^J (1 _ Q,2)s/2 



c+'i 00 



.T a2 0F 
'm2{2TTi) Jc-too 
( 3-s 4 



ds 



X 

Where we used that 
2^1 



-zmV c(^) r(^) 
2^ry' (i-a2)r(^) 



(60) 



l-a2 



1 



s 1 

2' 2' 2' 1-^2 



a2 - 1 



(61) 



In order to compute L\_i we can close the contour of 
integration by a half circle at Re(s) — *■ 00 since the con- 
vergence of the integrals is controlled by the ratio mjT 
which is taken to be less than one. The integral over this 
half circle vanishes. L\^2 are given by the residue of the 
poles of r[(3 — s)/2] and r[(2 — s)/2, which are located 
at s = 2fc + 3, and at s = 2(fc + 1), fc = 0, 1, 2 . . ., respec- 
tively. Choosing c > 1 we can avoid the pole of C(s) at 
s = 1. Working with this choice we get. 



2^\a^-l)l^^k\ 



2k+2 



! \2tiTJ 



X C(2fc + 3) 



2i^i(l,-fc,l/2-fc,^) 



L2 = 



I 

2^ 



2 



^a^ - 1 
X C[2(fc + 1)]. 



r(i/2-fc) 
1 



1/2 oa 



k=0 



1 



-( — ) 

)'= \2ttTJ 



2k+l 



(62) 



For L3 we notice that the contour of integration can be 
closed by a half cicle at Re(s) — > —00 since, as we will 
show, the contribution will be proportional to T/m which 
we take to be larger than one. The integral over this half 
circle vanishes. The integral is given by the residue of the 
poles of r[(s — 2)/2] which are located at s = —2k, k = 
0, 1, 2 . . .. However since C(— 2fc) vanishes for fc = 1, 2 . . ., 
the only pole that contributes is the one at s = 0. We 
thus get 



T 

I I — 

m 



\J a^ 



1 



(63) 



E- 



{2ttT} 



:C(^)- 



By changing a;„ —i^n in Eq. ([54|) . or equivalently, 
(59) ^ ~^ (ESI) and ([57)l . it is straightforward to 

show that for the modes with n < the contribution from 
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the the logarithmic term in Eq. (|57|) cancels. Therefore, 
the contribution from the modes with 71 7^ is just twice 
the above discussed contribution from the modes with 
n > 0. The remaining term to compute is the one coming 
from the mode with n = 0. It is easy to show that the 
delta function in this case does not have support and thus 
this contribution vanishes. 

Writing all together, the final result expressed as an 
explicit power series in the ratio m/2'iTT can be written 
as 



87r2 



j a 1^ \/q;2 — 1 — 



1 



1 



7rT\ 

TO J 
/ TO 

3 



C(2) 



f^rc(4) 



(a2 - 1) V27rT/ 



a" -1 
1 



'(55.)^ CP) 



2 V27rT 



+ 

+ 



"C(5) 



4(a2 - 1) V27rr 
(Sa" - lOa^ + 15) / to n 6 



16(a2 



2-kT 



C(7) 



(64) 



As we show in Appendix II, this result coincides with the 
one computed with the standard method. 



D. Infrared limit 

As can be seen from Eq. ([57)) [which is written before 
computing the integral over u for which we have assumed 
a > 1] in the infrared limit (a = 0), the self energy 
does not depend on a. In order to compare the result in 
Eq. ([5^ with Eq. (PS)) . which is computed in the infrared 
limit, we need to subtract from Eq. (|52|) the vacuum con- 
tribution, since the MST does not explicitly separate this 
from the thermal contribution. From Eqs. PH)) and (fTT|) 
the vacuum contribution to 11 in the infrared limit is 



= (0,p^0) 



8(27r) 



where we have extended the integral to d— dimensions. 
Notice that in this case, in order to make contact with 
Eq. (fTU]) . the limit d ^ \ will be eventually taken. Ex- 
plicit evaluation of Eq. (p5|) yields 



0) 



A2 



8(27r) 



-1e\. (66) 



Therefore, the thermal contribution is obtained by sub- 
tracting Eq. ([5^ from Eq. (|52p and this is given by 



n^(o,p^o) 



n(o,p-> 0) -n^^'^(o,p^ 0) 

A2 fTTT 
TO 



4(27r 



In f ™ ) 
\2TJ 



IE 



cm 



87r2r2 j 



(67) 



which coincides with Eq. (p5|l . 

At this point it is important to underline that the in- 
gredient making possible that the standard procedure re- 
produced in Sec. El] and the MST described in Sec. IV] 
lead to the same result is the implementation of the pe- 
riodicity of the expressions -by appealing to the fact 
that the external frequency is discrete- before the ana- 
lytical continuation to arbitrary complex values of the 
external frequency is taken. Also, when the Feynman 
parametrization is used and afterwards the periodicity 
implemented, the procedure leads to the well known re- 
sult, provided the integration domain for the Feynman 
parameter x is x € [0,1]- Nevertheless, as is discussed 
in Refs. P, Q, before the analytic continuation, the inte- 
grand is symmetric about x — 1/2 and thus it is seem- 
ingly possible to get the Feynman integral as twice the re- 
sult when X G [0,1 /2] . We show in Appendix III that this 
introduces the extra complication of a spurious end-point 
singularity and thus leads to the well known mishaps with 
the use of the Feynman parametrization in the ITF. 

We also point out that, as mentioned in Ref. [ij, for 
practical purposes, the result in Eq. (|67p . that is to say, in 
the infrared limit, can be directly obtained from Eq. ([1]) 
by setting poup = right from the start. In this case, 
in the context of the MST, the mixed frequency sum in 
Eq. collapses to the computation of the contribution 
of the n = Matsubara frequency and the definite sign 
frequency sums in Eq. (j27p can be condensed into a single 
sum over positive definite frequencies. For calculations 
involving propagators raised to higher powers, where one 
seeks an answer in the infrared limit, this simplification 
makes the MST to be a rather convenient technique, par- 
ticularly in the high temperature limit T ^ m since it 
gives the final answer in terms of a series in m/T. We 
proceed to show that this is the case when computing the 
self-energy of a scalar particle interacting with charged 
scalar particles in the presence of an external magnetic 
field. 



V. APPLICATION: SCALAR SELF-ENERGY IN 
A MAGNETIC FIELD 

In the Standard Model after symmetry breaking, there 
is an interaction term of the physical Higgs </> with the 
charged ones of the form Xip^ip(f>. In the presence of an 
external magnetic field, the propagators for the charged 
modes are affected, becoming, in the weak field limit and 
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at finite temperature 



2{eBfkl 



{eBf 



(68) 



where eB is the coupHng of the charged scalars to the ex- 
ternal magnetic field. One of the diagrams contributing 
to the physical Higgs self-energy at one-loop, depicted in 
Fig. [U is given explicitly by 



(fik 



(27r)3 

X D^{uo.^,k)D^{uo,,-uuk-p). (69) 
To lowest order in e_B, this self-energy becomes 

n^(^,,p) = X^tY, I ^ - (ei3)M/3i 

+ /i3-2A;i/4i -2(fc-p)i/i4]}, (70) 
where we define 

1 



(71) 



When interested in describing the infrared properties of 
this self-energy, we look at the infrared limit which, as 
previously discussed can be obtained in a straightforward 
manner by setting poi,P = 0. In doing so, we get 

n^(0,p-0) = X'tY: J ^{l20-2{eBf 

X [/4o - 2A:i/5o] } . (72) 



Notice that the functions /„o are all related through 



InO — 



i-iy 



gn 



(n - 1)! dim^y 



'10- 



(73) 



The MST technique discussed in Sec. |TV] can be gener- 
alized to the computation of Iio and from Eq. ([75)1 to 
all the expressions involved in Eq. ((72)) . The interested 
reader is referred to Refs. [1] for details and the result 
at high temperature and to lowest order in the magnetic 
field strength is 



n^(0,p->0) A' 



(47r)2 
T (eB) 



1 1 / M 



87rm 



64 



T 



1 C(5) 

7rm5 ^ 167r6 



(74) 



where we have not subtracted the vacuum contribution. 

We emphasize that the MST is suited to obtain an 
expression such as Eq. (I74p . namely, an expansion at 




FIG. 1: One-loop Feynman diagram contributing to the self- 
energy of the physical Higgs, represented by the single line, 
interacting with the charged Higgs components, represented 
by the double lines, in the presence of an external weak mag- 
netic field. 



high temperature starting from the original expression 
for the self-energy in the ITF. This is so since the sum 
over Matsubara frequencies and the integration over the 
spatial components of the momentum can be carried out 
together in a single step, in a straightforward manner, 
right from the very beginning. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have shown that the MST is a well de- 
fined method to compute Feynman integrals in the ITF 
of finite temperature field theory. The MST is particu- 
larly useful to find the explicit result as a series in a small 
parameter, for instance, the ratio m/T at high tempera- 
ture and in the infrared limit. The method calls for the 
use of the Feynman parametrization which in the past 
has been linked to problems in the ITF when the ana- 
lytical continuation from discrete Matsubara frequencies 
to arbitrary complex values is performed. We have also 
shown that these problems are not endemic to the Feyn- 
man parametrization and have traced back their origin to 
(a) not implementing the periodicity of the expressions 
before analytical continuation and to (b) changing the 
domain of integration in the Feynman parameter from 
X S [0, 1] to a; S [0, 1/2] which introduces a spurious end- 
point singularity. We have explicitly shown that when 
using the MST, and the calculation is properly carried 
out, it leads to the same result obtained by means of the 
standard technique in the ITF in the infrared limit. In 
particular, we have shown the need to take into account 
the correction term to the usual Feynman formula, in or- 
der to consider the case when the sum of denominators 
vanishes, and that this term is the source of the full a 
dependence of the result, in agreement with Ref. 0]. The 
usefulness of the method is illustrated by the computa- 
tion of the one-loop self energy in the standard model 
of the physical Higgs field interacting with the charged 



12 



components in the presence of a weak external magnetic 
field, in the infrared limit. 
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Appendix I: Evaluation of n(poi,p) without 
implementing the periodicity 

In this appendix, we aim at furthering the argument 
on the importance of having implemented the periodicity 
in the function coth in Eq. ([5]), to achieve the proper 
analytic continuation to arbitrary complex values of poi 
in the evaluation of n(po/,p)- We show here that when 
Poi is not taken initially as i times a discrete Matsubara 
frequency, then, when poi po, where po is a continuous 
arbitrary complex number, one is bound to obtain an 
spurious term which needs to be canceled precisely by 
the addition of the quantity tts of Ref. Q • 

Without using that poi takes on discrete integer values, 
instead of arriving at that equation we would have 



A2 f (fk 



E 



rs 



^8EkEk-pJ sEk-{rEk -p + Pol) 
sEk \ f rEk-p + poi ' 



coth 



2T 



coth 



2T 



(75) 



which can be rewritten as 



A f 

^{PohP) = -y / 

X E 



r.s— ±1 



(27r)3 
rs 

SEkEk- 



coth ( — - ) — coth 
2T 



sEk - [rEk-p +Pgi) 
rEk-p " 



2T 



csch' 



coth 



V 2T ) 



■ coth (ffr] 



(76) 



where we used coth(a + 6) = coth(a) — csch^(a)(coth(a) + 
coth(6))^^ to separate the dependence on Ek-p and pqi 
in the second hyperbolic function. Note that, compared 
to what we had in Eq. ([5]), we now have a third term as 
a result of not fully exploiting the periodic properties of 
the functions involved. 



We now concentrate in the last term in Eq. (|76|) and 
show that, according to Ref. Q and in the limit when 
p — > 0, this corresponds to minus the function one needs 
to add to correct the result. Let us then call lixiPvihP) 
the contribution from the aforementioned term 



nx(po;,p) 



A2 
2 



d^k 
(27r)3 



E 

r,s— ±1 



SEkEk- 



sEk - (r Ek-p +Poi) 
csch^ 



coth(^%^) +coth(£aL) 
where, upon summing over s we have 

A2 r d^k 



(77) 



UxipohP) = -- 



E 



1 



{rEk-p +PQiy 
csch^ 



El 



2T 



coth 



2T 



coth(H^) 



(78) 



In order to integrate out the angular contribution, we can 
perform the momentum shift k — p ^ k so that all the 
angular dependence will be in the coefficient rather than 
in the hyperbolic functions. This allows for a straight- 
forward integration and we arrive at 



^x{poi,p) 



A^ 
2 



dk x ^ 

(2^ ^ 

CSCh^ (^) 



— rfc 

1e7 



coth 



coth (ff.) 



rEk 
2T 

1 , / pIi - P^ + '2rEkPoi - 2kp 
2kp i ^2 _ „2 



Poi - + '^rEkPoi +2kp 



(79) 



We now proceed as in the main body of the paper af- 
ter Eq. (O. We take the analytical continuation in poi 
from discrete imaginary values to arbitrary complex ones 
Pol ^ Po- Since the result depends on how the limit is 
explored, we first set po = oip. To analyze the behavior 
near the origin, we expand the function coth and the log- 
arithm around p — Q and, up to linear terms, we obtain 



Xlx{ap,p) 



p^O 



A^ 
2 



dk 



Ef 



-rk \ 



(2^f L.^\SEkp) 





(rEk\ 




y2T ) 



In 1 - 



2k 



k + raEk 



(80) 
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which, after summing over r gives 



Uxiap,p) -"^ 



In 



2 JO 
k 



dk 
aEkJ 



Yk 



esch 



2T 



(81) 



Now, just as we did in Sec. Illli we are interested in 
having an expHcit functional dependence on a of IIx- 
We can then easily extract the thermal contributions 
thereby knowing how Wx modifies 11, as was discussed 
in Eq. (I12p . For this purpose, it is convenient to note 
that the term in the square brackets of Eq. ([8T|) can be 
written in terms of a partial derivative {dkCotlvEk = 
kE^^ csdr? Ek) , so that we can complete a total derivative 
through integration by parts, to have 



„ . . P^o A2 dk c a 



A2 
2 

d_ 

dk 



coth 



Ek 
2T 



In 



k — aEk \ ' 



2Ek{k-^ -a-^Ei) 



k + aEh J 
/ I 

coth 



Ek 
2T 



(82) 



Finally we can separate the vacuum and the thermal con- 
tributions using the identity in Eq. (jlip , so that the ther- 
mal part is 



nx(ap,p) 



A2 



2(27r)2 



dk 



n{Ek 



Ek [k^-a^ElY 



(83) 



The function li^ in Eq. (|83|) is precisely limp^o t^s {ctp, p) 
found in Eq.(30) of Ref. [3|, but with the opposite sign. 
We can see that in the event of not implementing the peri- 
odicity, as we have analyzed in this appendix, inevitably 
we will end up with a contribution stemming from the 
extra term IIx. The situation is corrected, as noted in 
Ref. [1], if one adds a function that behaves just as XI^ 
in the limit considered. This turns out to be an impor- 
tant observation, since we are presenting evidence that 
neglecting the implementation of the periodicity in the 
external frequency is linked to the need of such correct- 
ing function. Further developments on this argument are 
presented in the rest of this work. 



Appendix II: Evaluation of n^(ap,p) for p — > and 
arbitrary a 



We start from Eq. p2|) rewriting it as 



n iap,p) 



2(27r)2 J, 



Ek El - a'^m^ 



(84) 



where a'^ = j^^- We follow again Ref.p^ and use the 
identity in Eq. ([15]) into Eq. ([M]), that is 



Il^{ap,p) 



X' //-'^ /■ ... 1 



d'^k- 



1 

2(27r)2 iiT J ^ '^El~a'^m^ \ 2£^ 
_ 1 1 



(Ek)^ + (27rnT)2 



(85) 



where we have written the integral in c?-dimensions. The 
first structure in Eq. ([55]) is 



Ji = -M' 



3~d 



X' 1 



d'^k 



1 1 



2(27r)2 47r J 2Ek El - a''^m^ ' 



(86) 



Carrying our the angular integration and upon the 
change of variable z — , we get 

Jl = _^3-d^'i-3--^-Ll-^ / dz(l_z)4-l 



(87) 



2(2^)2 87rr(f)7o 

Using the identity 

2F,{a,b,c;z) = ^^J}'''^ C dt{\ - tf-^-^ 



r(5)r(c-5) Jo 
X t^-^{l- zt)-", 

where 2F1 is the hypergeometric function, we get 



Jl = -/i^""™ 



A2 r(3^) 



2(27r)2 8^ r(|) 



X 2F1 ( l,^-,-;a 



(89) 



For the second structure in Eq. (|85p. a similar proceedure 
leads to 



J2 = 



3-d 



-^^r(2-^ 

2(27r)2 47r V 2 



, d Lu?, + a'^m'^ 

X l,2--,2; " ^ 

2 w„ -I- m-^ 



(90) 



Note that for the term n = in Eq. ((90| . the argument 
of the hypergeometric function becomes independent of 
m and T. The result for Il{ap,p) is thus 



n^(ap,p) = Ji + J2. 



(91) 



In order to veryfy this result in the limit a = {a' = 1) 
we recall the identity 



2^1(0, 6, c; 1) 



r{c)r{c - b - a) 

r{c~b)r{c~ a)' 



(92) 
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that can be used to write 



.^0 1 /M^3-..(^-l)/^^(^)^(f-l) 

^ 2(27r)2 Stt Vm/ 

.=0 . . 1 ^^/^r(2-i)r(i-i) 

J2 ^ M 



r(|) 



2(27r)2 47r r(f) 

OO 



(93) 



n— — ctt 



Using the procedure as in Eqs. (|20|) - ((23)) to obtain the 
high temperature Umit, we get for J 2 

^ A2 1 3_, ^'^/2r(2-f)r(f-i) 

2(27r)2 4^^ r(f) 
X ( Tm'^-^ + 2T{2TTT f-*C{4: - d) 

X +2T - 2] {27rTf-'^m^C{6 - d)\ (94) 



Taking d ^ 3 — 2e and a = 0, the result in the infrared 
hniit is 



n^(o,p^o) = 



4(27r) 
m2C(3) 
87r2T2 



h hi — 

m V2r/ 



7s 



(95) 



which coincides with Eq. 

We can also use the former analysis to give an explicit 
expression for the a dependence of the self-energy in the 
high temperature limit. We first separate from Eq. ((84|) 
all a dependence. In terms of the parameter a', we get 



n(ap,p) 



A2 



2(27r)- JO 

Ho + n„, 



dk 



njuJk) 



(a'2 - l)m2 



tut — a 



(96) 



from where we get 

A2(«'2_i)ir(i) 



Ji 



J2 



2(27r)2 4r(|) 
A2(a'2-l)l sin-i(a') 



3^2i^i(l,l,2;«") 



2(27r)2 2a'VT^^ 



X^a'^ ~ I) nTm- 
2(27r)2 2 



2^^1(1, ^,2; a'2) 



+ 2j2xlixl + l)-^2F^{l,-,2;^-^ 

71=1 " 

= Jl+J2. 



(99) 



where x„ = m/2TrnT. Notice that Ji in Eq. ([M]) yields a 
T-independent term and therefore contributes only to the 
vacuum part. This can be shown to correspond to con- 
sidering the pole of C{s) at s = 1 in L2 given in Eq. ([50)1 . 
We can thus ignore this term. 

In the high temperature limit the parameter a;„ <IC 1, 
thus we can perform a series expansion in J2, yielding 



J2 



A2 



ttT 



1 - a 



+ 



a'- 



o? - 1 
1 



(a2 - 1) V27rT/ 



rrw (2^) C(6) 



(a2 - 1) 
a2 - 1 



1 / m \2 

2 (2^) 



4(a2-l) 12^^ J ^^^^ 
(Sa'^ - 10a2 15) / m 
16(a2 _ 1)2 



C(7) 



(100) 



which coincides with Eq. 



Appendix III: in the ITF 



where 



n(a:p,p)c 



A2(a'2-l)m2 r^^niu^ 



2(27r) 



9 /9 9 ■ 



(97) 



Notice that the above integral can be obtained from 
Eq. ([M|) taking d = 1 — 2e and with the changes 



1 1 

4^ ^ 2 



(98) 



In this appendix we show that the function II^ found 
in Ref. 0] emerges in the ITF making use of the Feynman 
parametrization only when the limits of integration are 
replaced from x G [0, 1] to cc G [0, 1/2]. 

,f.l3 



We start from Eq. (3.33) in Ref. 

A2 r d^k 



(27r)3 



dx 



d 



dim? 



E 

r=±l 



coth ^(r xpoi + J/2 ) 



where y is defined as 
y = El + x{Ef 



El) 



x{l 



x)pIi. 



(101) 



(102) 
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It is worth noticing that in Ref. [s*] the change of variable 
k — xp ^ k is performed in Eq. (jlOip . but this change 
is not allowed in this case since the integral is divergent, 
unless the divergence is regulated by using for instance, 
dimensional regularization. 

Using the identity in Eq. (|48]) into Eq. (llOip . we get 



r— ±1 



(2^ 



dx 



d_ 

dx 



coth|(r xpm +2/2) 
22/(rpoi + #) 



The integral over x becomes trivial and when evaluating 
in the integration limits a; = La; = we obtain Eq. 
This is what is done in Ref. [l|] which leads to the cor- 
rect result, provided the periodicity in coth is imposed, 
as discussed in Sec. IIIIl However, if we instead follow 
Ref. 3 and use that the integrand is symmetric about 
a; = 1/2 and thus that the integral over x in the interval 
X E [0, 1] is twice the integral in the interval x € [0, 1/2], 
we get 



A2 

' 4 



E 

r=±l 



d^k 
(2^ 



dx 



dx 



coth I (r xpoi +2/2) 



(104) 



Notice that Eq. (|104p is valid when poi is imaginary and 
discrete, since only in this case, coth is periodic. Evalu- 
ating the integral over x in Eq. (|104p we get 



^xiPOhP) 



X2 



r=±l 



d^k 
(2^ 

coth ^Ek 



+ 



Ek{2Tp^iEk + Ei_^^Ei-pi,) 
coth ^|(r xpoi +yi 



2y{r poi + ^) 



x=l/2 



(105) 



where the first term results from evaluating in the lower 



limit of the 2;-integral and in the second one we have left 
indicated that x is evaluated in 1/2. Notice that when 
completing the square in the denominator of the first 
term in Eq. (jlOSp . this becomes identical to the result in 
Eq. ([6]), which is the correct result, thus leaving Eq. (|105[) 
with an extra term, which in fact, as we proceed to show, 
corresponds to the function —Us in Ref. To show this 
we must carry out the angular integration in Eq. pOSp . 
Defining 



(103) IlAipoi,p) 



X2 

4e 



r=±l 



d^k 
(2^ 



coth I (r xpoi + J/2) 



[2ryhoi + El_^ - El - (1 - 2x)pl\ ' 

(106) 

where x should be evaluated in 1/2. Upon the change of 
variable k — a:p k, the dependence of the angle inside 
the function coth is removed and we get 



nA(poj,p) 



A2 ^ /■ d^k 
T ^ y (2^ 

coth I (r xpQi -f 05) 
(l)^2r(f>ipoi - 2k • p] 



(107) 



where (f) = k'^ + — a;(l — a;)(pgj — p^). The remaining 
angular integration is readily performed and the result is 



nA(po;,p) 



4(27r) 



kdk 



Poi(, 



kp 



n{xpoi -I- 02) 



. In 

ip \^po/02 - kp 
' n{-xpoi + (j) 



x=l/2 

(108) 



where n is the Bose-Einstcin distribution. Notice that if 
in Eq. (llOSp we use that po; is purely imaginary and dis- 
crete, the function Ha vanishes. However, if po/ is analyt- 
ically continued to arbitrary complex values, the correct 
result is obtained only by the addition of the function H5 
found in Ref. ji], which exactly cancels Ha- 



[1] P. S. Gibosky and B. R. Holstein, Z. Phys. C 47, 205 
(1990). 

[2] P. F. Bedaque and A. Das, Phys. Rev. D 45, 2906 (1992). 

[3] H. A. Weldon, Phys. Rev. D 47, 594 (1993). 

[4] A. Sanchez, A. Ayala and G. Piccinelh, Phys. Rev. D 75, 

043004 (2007); A. Ayala, A. Sanchez, G. Piccinelh and 

S. Sahu, Phys. Rev. D 71, 023004 (2005). 
[5] D. J. Bedingham, Dimensional regularization and 

Mellin summation in high temperature calculations, 



|hep-ph/ 0011012| 

[6] G. Metikas Uniqueness of the thermal effective potential, 

'hep-th/9910063 
[7] T. S. Evans, Z. Phys. C 41, 333 (1988). 
[8] T. S. Evans, Derivative expansions of Euclidean thermal 

effective actions, hep-ph/9808382 
[9] H. A . Weldon, Phys. Rev. D 28, 2007 (1983). 
[10] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 
[11] Aa Kilbas, Handbook of Mellin Transforms, (Cambridge 



16 

Scientific Pubfishors, 2006). Company, 1995). 

[12] M. E. Peskin and D.V. Schroeder, An Introduction 
to Quantum Field Theory, (Addison- Wesley Publishing 



